
function z = mvn_rnd(mu,sigma)
z = mu + chol(sigma)'*randn(length(mu),1);